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ABSTRACT 

Direct imaging of the HR8799 system was a major achievement in the study of 
exoplanets. HR8799 is a 7Doradus variable and asteroseismology can provide an 
independent constraint on the inclination. Using 650 high signal-to-noise, high 
resolution, full visual wavelength spectroscopic observations obtained over two 
weeks at Observatoire de Haute Provence (OHP) with the SOPHIE spectrograph 
we find that the main frequency in the radial velocity data is 1.9875 d -1 . This 
frequency corresponds to the main frequency as found in previous photometric 
observations. Using the FAMIAS software to identify the pulsation modes, we 
find this frequency is a prograde £—1 sectoral mode and obtain the constraint 
that inclination i >40°. 

Subject headings: stars: oscillations — stars: variables: other — stars: individual 
(HR8799) 
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Introduction 



The imaging discovery of the three (IMarois et al. 1120081 ). and now four fjMarois et al. 



20101 ) planets around HR 8799 is a significant achievement in the search for and study of 
planets orbiting other stars. For the first time, the thermal emission of planets in orbit 
around another star has been unambiguously detected. 

The dynamical evolution of a planetary system is complex. From the basic planet forma- 
tion assumption that planets form by the core accretion or disk instability scenario in a disk 
along the star's equatorial plane, systems can suffer drastic changes; planet-planet perturba- 
tions, interactions with a disk or stellar encounters can cha nge a planet's orbital inclination 



its se mi-major axis, its orbital eccentricity or even eject it ( [Raymond. Armitage fc Gorelick 
2010). In the case of HR 8799, it is not impossible that a close encounter occurred prior 



to planet formation, since its relatively high galactic velocity compared t o the Columba as- 



sociation and its far distance away from the other association members (iHinz et al. 1120101 ) 
are suggesting that it may have been kicked out and is probably moving quickly away from 
its birth place. Such an encounter could have tilted a disk relative to the star's equatorial 
plane and induce perturbations that may have led to planet formation where planets would 
have a non-negligible orbital inclination relative to the star. Radial velocity searches have 
confirmed such chaotic behaviors by detecting system s where planets are orbiting well away 
from the star's equatorial plane ( jTriaud et al. 1 12010| and references therein), although for 
close -in extrasolar planets this misalignment could also be caused by the Kozai mech anism 
(e.g. IWu. Murray fc Ramsahai 1 120071 : iFabrvckv fc TremainelboOTl : Iwinn et al. IlioOa ). For 



the wide HR8799 planets the Kozai mechanism is not operational (even if the system had 
a stellar companion in the past); finding a misalignment between the star and the planet's 
orbital plane would be a sign of a significant dynamical interaction in the system past. 



Marois et al. I ( 120081 ) have suggested that the HR 8799 planets are in a similar orbital 
plane with a low inclination and have mostly circular orbits. This is because the detected 
orbital motions are close to their expected face-on circular orbit values, the orbital motion is 



pole- 


on). Dynamical analyses (e.g. 


Reidemeister et al. 


2009 




Fa 


jrvcky & Murray- Clav 


2010 




Moro-Martm. Rieke & Su 


20ld 


Marshall. Horner & Carter 


2010) have confirmed 



that the planets are mostly in the same plane with small eccentricities, although such fits 
are still very uncertain due to the limited amount of orbital coverage available. In addi- 
tion, the planets also all orbit in the same counter-clockwise orientation, further support- 



1 Based on observations obtained at the Observatoire de Haute-Provence which is operated by the Institut 
National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France. 
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ing the idea that the y formed in a disk , simi lar to the Solar system planets. Assuming 
a circular orbit for b, lLafreniere et al. I ( 120091 ) found % ~13°-23°; while attempting a co- 
here nt analysis of various porti ons of observational data on known components of the sys- 
tem, iReidemeister et al. I ( 120091 ) concluded that % should range between 20° and 30°. Also, 
Spitzer observations of HR 8799's co mplex debris disk suggest that any inclination angle 
larger than ~25° sho uld be excluded (ISu et al. 1120091) . Using a statistical distribution of 
star's rotation speed (Rover. Zorec fc Gomez 1120071 ). HR 8799 with its 37.5±2 kms" 1 Vsinz 
( iKaye fc Strassmeier 1 ll998h would be consistent with an inclination of ~ 23.5° if it is an A5 
star or ~ 18.5° if it is an F0 star (HR 8799 spectral classification is uncertain mainly due 
to its low metallicity that is affecting it's broad band colors). Such a determination is of 
course statistical and a direct star's inclination determination is required for a meaningful 
comparison with the estimated planet orbital plane inclination. 



HR 8799 is an intrins ic photometric and spectroscopic variable ( jRodriguez fc Zerbi 



1995 




Mathias et al. 


dZerbi et al. 


1999 





. 120041 ) . It has also been confirmed as a 7Doradus (Dor) variable 
jZerbi et al. 19991 ). The 7 Dor stars are late A to early F stars whose pulsat ions are driven 



by a flux-blocking mechanism at the base of their convective envelope (e.g. iDupret et al. 



20041 ). The 7 Dor nature of HR 8799 offers a unique opportunity to estimate its inclination 
via an asteroseismic analysis of the observed q-m odes. In a previous asteroseismic analysis 
using photometric frequencies iMoya et al. I ( 120101 ) has shown that an age determination for 
the system, which allows to discriminate between planets and brown dwarfs, is a difficult task 
with the current information and discussed the importance of an inclination determination 
since the equatorial rotational velocity can be used in constraining the age of the system. 
In this letter, we present a spectroscopic asteroseismic analysis and also obtain limits for 
HR 8799's stellar inclination. The data have been acquired from an extensive multi-site 
ground-based high-resolution spectroscopy campaign. Section 2 describes that data used in 
this letter. Section 3 discusses the pulsation mode identification and the determination of 
the stellar inclination and Section 4 describes the conclusions. 



2. Observations 

In this letter we examine data collected from Observatoire de Haute Provence (OHP) 
using the SOPHIE spectrograph. We have 650 SOPHIE observations taken from 28/09 - 
12/10/2009. They cover a wavelength range between 3875-6940A with a spectral resolution 
of R~75 000. These data are the largest single-site subset of more than 2000 observations 
taken during an intensive spectroscopic multi-site campaign devoted to this star. Only the 
OHP spectra will be examined in this publication as study of the complete combined data 
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set of observations will be left until the detailed publication. 



The data were treated using the local reduction packages written specifically for the 
SOPHIE spectrograph data. This produced reduced, merged and automatically continuum 
normalized 1-D spectra which were further manually normalized using a synthetic spectrum 
for guidance. 

Instead of using single spectral lines, we examined the cross-correlation profile (CCP). 
This profile is obtained using the normalised 1-D spectra cr oss-correlated with a lin e mask. 
The line mask was produced with the SYNSPEC software (IHubeny &: Lanz 1119951 ) using a 
Te//=7500 K, log(g) = 4.0 and Z=-1.0. All the telluric and H/He affected areas were removed 
from the spectral line for which we estimated the equivalent widths. These equivalent widths 
were used for the weightings given to the delta-functions in the cross- correlation. An analysis 
of the variation of several isolated spectral lines gave similar results to those outlined for the 
CCPs below. Since the signal-to-noise ratio of the individual spectral lines is much lower 
than that of the CCPs only the CCPs are discussed here. 

The resulting line variation was examined using the first moment of the line profile 
which is similar to the radial velocity of the line (shown in Figure 1: left). A Fourier analysis 
of the first moment gives two strong frequencies 1.9875 d _1 and 1.7475 d -1 (shown in Figure 
1: right). These frequencies by themselves produce a good fit to the observed r adial velocity 



variat ions. They closely correspo nd with those found in photometry by both IZerbi et al. 
(|l999h and lCuvpers et al. I J2006h . 



3. Pulsation mode identification 

It is possible to identify further frequencies present in the data but upon examination 
it is observed that extraction of any combination of the other frequencies present does not 
alter the results for the 1.98 d^ 1 frequency. The focus for this letter is on the inclination 
restrictions possible from the 1.98 d _1 frequency, hence we will continue without examining 
any further frequency information, and leave that for the detailed publication based on the 
whole dataset of the campaign. 

To determine t he degree £ and azimuthal order m of the pulsation mode using the 
FAMIAS softwar^l ( Zima 2008 ) we extract the Fourier parameters, i.e. the zero point, 



the amplitude and the phase distributions across the CCPs, for the two frequencies. The 



2 The software package FAMIAS developed in the framework of the FP6 European Coordination Action 
HELAS (http://www.helas-eu.org/! 
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Fig. 1— OHP Radial velocity data and 1.9875 d" 1 and 1.7475 d" 1 fit for the first half of 
the OHP data (left). Frequency analysis of the radial velocity of the OHP data (right) with 
lines marking the frequencies determined. 



Fourier parameters of 1.9875 d -1 are shown as a solid line in the panels of Figure 2. The 
Fourier parameters of 1.7475 d _1 were not sufficiently well matched by the FAMIAS software 
to consider this pulsation mode identified or constrained. This is probably a result of other 
pulsation frequencies present in the data that have not been detected or removed, and that 
have sufficient amplitude to affect the parameters of the 1.7475 d -1 mode. 

The FAMIAS software uses a first order Coriolis force approximation, limiting the pul- 
sation models that can be used to fit the Fourier parameters to those respecting v < 1, where 
v is the so-called "spin-parameter". This parameter is defined as v = 2 fn/ fcorot, where 
is the rotation frequency and f coro t is the pulsation frequency in the corotating frame. f coro t 
is connected to the observed pulsation frequency f Q b s through: f corot = fobs - m fn and, hence, 
changes for different values of m of the model being fitted. As for fo, its changes depend 
on the inclination used by the model. On the other hand, one could wonder how significant 
is the effect of stellar rotation w hen v > 1 and if the fit of the Fourier parameters could 



Corio lis force on high order 



n) 
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2003 


) suffer 



reduced changes from increasing Coriolis force. Using this study we can estimate more ap- 
propriate v limitations for the FAMIAS models of the 1.9875 d -1 frequency. For the values 
of m=-3:3, limits of v <8 for m=3, v <5 for m=2, v <2 for m=l and v <1 for m <0 
are reasonable. Usin g these v limitations and the estimated stellar parameters R*=1.5R 
( Gray fc Kaye 1999h . Vsini=39.5 kms" 1 ( i.e. a value within the uncertainties given by 



this study and by iKaye fc Strassmeier l ll998l) we can place restrictions on the inclination for 
which the FAMIAS models are valid for each value of m, they are i= 60°, 37°, 32°, 32°, 15 c 
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0° and 0° for m— 3, 2, 1, 0, -1, -2 and -3, respectively. 

These model limits imply that we cannot test all the inclinations between 0° and 90°. 
However, one has to consider that an inclination lower than 5° is physically impossible, since 
it would imply an equatorial rotat ion velocity higher t han the break-up velocity when using 
an estimated mass of M*=1.5 M (IGray fc Kaye Ill999l ). Moreover, models with m=l, 2 and 
3 are limited to i >16°, 33° and 54° respectively to be physically possible (f CO rot > 0). A 
summary of the limitations on the inclination of the pulsational axis is given in Table 1. 
From this table it can be seen which are the physically possible values of m and inclination 
that we are unable to test due to the FAMIAS model restrictions. However, these "holes" can 
be mitigated by the fact that at lower inclinations than those imposed by the model limits, 
the v value, and hence the Coriolis force, becomes large en ough that the pu lsations surface 
deformations begin to be limited to an equatorial waveguide (ITownsend l l2003h . In these cases 
just the regions about the equator are varying, but we would be viewing the star from low 
inclinations. This means that only a small area of the visible surface would be experiencing 
pulsation and the radial velocity amplitude would be expected to be low when compared with 
observed radial velocity amplitudes in other 7 Dor stars. To give an example of this, if i=20° 
and the mode was m=0 then approximately one quarter of the full surface variability is visible 
since more than half the pulsation amplitude is constrained to within 35° of the equator. 
In contrast, the 1.98 d -1 frequency in question has an amplitude of 1.09 kms -1 which is 
comparable to the radial velocity amplitudes of the strongest mode in other 7 Dor stars e.g. 
1.3 k ms -1 in 7Doradus (IBalona et al. Ill996l ). 0.35 kms" 1 in HD 49434 (lUvtterhoeven et al 



20081 ). 1.45 kms" 1 in HD189631 and 0.49 kms" 1 in HD40745 flMaisonneuve et al. 1 120101 ) . 



Therefore, it is unlikely that we are dealing with a strong Coriolis force and hence a high 
v value for a given m value. Because of this the inclination "holes" mentioned previously, 
which are associated with high v values, can be considered as unlikely solutions. 



Table 1: Table of constraints on the inclination from both the model limitations and physical 
constraints 

m value 3 2 1 0-1-2-3 

model i limits (°) >60 >37 >32 >32 >15 >0 >F 
physical i limits (°) >54 >33 >16 >5 >5 >5 >5 



All £ and m combinations from £=0 to £=3 were tested keeping in mind the restrictions 
from Table 1. Higher values of £ were not considered b e cause the 1.9875 d _1 is a strong 
photometric frequency found in the data of IZerbi et al. I ( 119991 ) and modes with £ > 3 are 



hard to detect in such ground-based photometric studies due to geometric cancellation effects. 
The input parameters to FAMIAS were permitted to vary as shown in Table 2. The two 
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best fitting modes with different £ and m values obtained with FAMIAS are shown in Figure 
2. The dashed-dotted line corresponds to an £=1 m—1 mode with x^ ed =30.0 and the dotted 
line to an £=2 m—-2 mode with %^ ed =43.9. This figure confirms that the best fitting mode is 
the £=1 sectoral mode as the £=2 m—-2 solution does not match the shape of the amplitude 
distribution well. The extremely small uncertainties in the zero point distribution cause 
any deviation from a perfect fit to rapidly increase the xled v& hie. Hence, the value of 
Xr e d=30.0 should be considered as a very good fit. On examination of the amplitude and 
phase distributions across the profile for all of the modes tested, only the £=l m—1 mode 
demonstrates a good match in shape for both and is considered conclusively the best solution 
for this pulsation frequency. The values for the input parameters to FAMIAS resulting from 
this mode are given in the last column of Table 2. The errors are estimated based on an 
examination of the Xred results. 

Table 2: Table of parameters for genetic algorithm search. Both the allowed range and best 
fit results, an £=1 m—1 mode, are shown. Uncertainties are estimated based on x 2 red results 



Parameters 


Allowed 


ranges 


Best fit 




min : 


max 




M* (M ) 


1.2 : 


1.8 


1.5±0.3 


R* (R©) 


1.2 : 


1.8 


1.5±0.3 


Inclination (°) 


5 : 


90 


65±25 


Vsini (kms -1 ) 


36 : 


46 


39.8±0.4 


intrinsic width (kms -1 ) 


5 : 


15 


8.3±0.5 


pulsation amplitude (kins" 1 ) 


0.1 : 


10 


0.659±0.05 


pulsation phase (2%) 


: 


1 


0.95±0.03 



Figure 3 shows the inclination dependence of the genetic algorithm search for the best 
fitting £—1 m=l mode. Based on Table 1 we can confidently use the FAMIAS results for 
i >32°. The best fit inclination has z=65°. For i= 35°, 40° and 90°, we find xl e d= 37 - 2 > 32 - 9 
and 31.07, respectively. It is obvious from Figure 3 that the Xred does not change much in 
the range i= 40° to 90° but it rises quickly at lower inclinations. From this we obtain our 
result % >40°. 

The FAMIAS software assumes that pulsation axis is aligned with the rotational axis, 
but this is not necessarily the case in all stars. However, such misalignment is unlikely in the 
case of HR 8799 because, in that case, some modulation of the pulsation properties by the 
rotational period would have been observed. In addition, it would have led to more complex 
amplitude and phase across the profile, and the fits that we have achieved would not have 
been as good. 




Fig. 2. — Zero point profile (top), amplitude distribution (middle) and phase distribution 
(bottom) for the frequency 1.98 d^ 1 is the solid line. The best fit £=1 m—1 mode is the 
dashed-dotted line with a xled °f 30.0. The next best fitting mode of those tested is the £=2 
m—-2 mode shown as a dotted line with x^ ed =43.9. 
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Fig. 3. — The inclination dependence of the genetic algorithm search for the best fitting £=1 
m—1 mode. The inset figure is a zoom in on the region of low Xred- For inclinations i>40° 
the differences in xled ls small. 



It is worth noting that no planetary reflex velocity is identified. This was expected 
since the amplitudes of the quite distant and low orbital inclination exo-planets would be 
only fractions of a meter per second which would be masked by the much larger (kms -1 ) 
pulsational velocities and is also beyond the precision attainable by our observations. 
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Conclusions 



We conclude that the stellar rotational inclination axis has a value i>40° based on 
identification of the 1.98 d -1 frequency as a n 1=1 m=l mode. Th is is the strongest pulsa- 
tion in both photometry (IZerbi et al. I Il999l ; ICuypers et al. 1 120061 ) and spectroscopic radial 
velocities. Through dynamical analyses it is suspected that the planets are mostly in the 
same plane with small eccentricities and that the p lanets orbit inclination axis is ~20°±10° 



( jReidemeister et al. 112009c lLafreniere et al. 1 120091 ). The current data suggests a misalign 



ment of Ai >20° between the stellar rotational inclination and planetary orbit axes, though 
more detailed pulsational analyses and better orbital fits are needed before this can be con- 
firmed. 



We thank A. Moya for helping discussion and useful advices. Wright acknowledges 
support from the Belgian Federal Science Policy (project MO/33/021). ANC acknowledges 
support from Comite Mixto ESO-GOBIERNO DE CHILE and from BASAL/ FONDAP 
project. PDC acknowledges financial support from the Fund for Scientific Research - Flanders 
(FWO; project G. 0332. 06). Portions of this work were performed under the auspices of the 
U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract 
DE- AC52-07NA27344. 

Facilities: OHP-1.93m (SOPHIE) 



REFERENCES 

Balona, L. A., Bdhm, T., Foing, B. H. et al. 1996, MNRAS 281, 1315 
Cuypers, J., Aerts, C, De Cat, P. et al. 2006, A&A, 499, 967 

Dupret, M.-A., Thoul, A., Scuflaire, R., Daszyska-Daszkiewicz, J., Aerts, C, Bourge, P.-O., 
Waelkens, C. & Noels, A. 2004, A&A 414, L17 

Fabrycky, D. C. & Murray-Clay, R. A. 2010, ApJ, 710, 1408 

Fabrycky, D. & Tremaine, S. 2007, ApJ, 669, 1298 

Gozdziewski, K. & Migaszewski, C. 2009, MNRAS, 397L, 16 

Gray R. O. & Kaye A. B. 1999, AJ, 118, 2993 

Hinz, P. M., Rodigas, T. J., Kenworthy, M. A., Sivanandam, S., Heinze, A. N., Mamajek, 
E. E. & Meyer, M. R. 2010, ApJ, 716, 417 



- 11 - 



Hubeny, L, Lanz, T. 1995 ApJ, 439, 875 

Kaye, A. B., Strassmeier, K. G. 1998 MNRAS, 294L, 35 

Lafreniere, D., Marois, C, Doyon, R. & Barman, T. 2009, ApJ, 694, L148 

Maisonneuve, F., Pollard, K. R., Cottrell, P. L. et al. 2010, MNRAS 

Mathias, P., Le Contel, J.-M., Chapellier, E. et al. 2004, A&A 417, 189 

Marois, C, Macintosh, B., Barman, T. et al. 2008, Sci 322, 1348 

Marois, C; Zuckerman, B.; Konopacky, Q. M.; Macintosh, B.; Barman, T.; 2010 arXiv 
1011.4918 

Marshall, J., Horner, J. & Carter, A. 2010, IJAsB, 9, 259 
Moro-Martm, A., Rieke, G. H. & Su, K. Y. L. 2010, ApJ, 721, L199 

Moya, A., Amado, P. J.; Barrado, D., Garcia Hernandez, A., Aberasturi, M., Montesinos, 
B. & Aceituno, F. 2010, MNRAS, 405L, 81 

Raymond, S. N., Armitage, P. J. & Gorelick, N. 2010, ApJ,711, 772 

Reidemeister, M., Krivov, A. V., Schmidt, T. O. B., Fiedler, S., Miiller, S., Lohne, T. & 
Neuhauser, R. 2009, A&A, 503, 247 

Rein, H., Papaloizou, J. C. B. & Kley, W. 2010, A&A, 510, 4 

Rodriguez, E. & Zerbi, F. M. 1995, IBVS, 4170, 1 

Royer, F., Zorec, J. & Gomez, A. E. 2007, 463, 671 

Su, K. Y. L., Rieke, G. H., Stapelfeldt, K. R. et al. 2009, ApJ, 705, 314 

Triaud, A. H. M. J., Collier Cameron, A., Queloz, D. et al. 2010, A&A 524A, 25 

Townsend, R. H. D., 2003, MNRAS 343, 125 

Uytterhoeven, K., Mathias, P., Poretti, E. et al. 2008, A&A, 489, 1213 

Winn, J. N., Johnson, J. A., Albrecht, S., Howard, A. W., Marcy, G. W., Crossfield, I. J. & 
Holman, M. J. 2009, ApJ, 703L, 99 

Wu, Y., Murray, N. W., Ramsahai, J. M. 2007, ApJ, 670, 820 



Zerbi, F. M, Rodrigez, E., Garrido, R., et al. 1999, MNRAS 303, 275 
Zima, W. 2008, CoAst, 155 



This preprint was prepared with the A AS IATgX macros v5.2. 



